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Abstract —As penetration of wind power generation increases, 
system operators must account for its stochastic nature in a 
reliable and cost-efficient manner. These conflicting objectives can 
be traded-off by accounting for the variability and uncertainty of 
wind power generation. This letter presents a new methodology 
to estimate uncertainty sets for parameters of probability distri¬ 
butions that capture wind generation uncertainty and variability. 

Index Terms —Wind power uncertainty, wind power variability, 
power system operations. 

1. Introduction 

W IND power generation (WPG) introduces variability 
and uncertainty in operational planning. As defined 
in m, variability of WPG is the random fluctuation of wind 
speed caused by physical processes in the atmosphere while 
uncertainty of WPG results from wind forecast errors. To ac¬ 
count for variability and uncertainty, approaches to robust unit 
commitment (RUG) la), chance constrained optimal power 
flow (CC-OPF) lO and distributionally robust CC-OPF (RCC- 
OPF) l[3|,p4l have been formulated and tested. However, the 
performance of these models depend on the accuracy of pa¬ 
rameters of probability distributions that define the uncertainty 
sets used; therefore, the uncertainty in distribution parameters 
must be accounted for. This letter makes two contributions: 
i) we use a data-driven analysis to relate intra-hour wind 
speed variability to the hourly-average wind speed and ii) 
using this relationship we construct uncertainty sets that can 
be interpreted in terms of variability and uncertainty of WPG. 

Bienstock et al 10 show that the CC-OPF formulation 
with Gaussian-distributed deviations of WPG and precisely 
known mean and variance can be extended to the RCC-OPF 
formulation, where the parameters of the Gaussian deviations 
(both mean and variance) fall within uncertainty sets. This 
RCC-OPF is implemented and tested on a large-scale system 
in n. In this letter, we interpret the set for the mean in terms 
of the wind power uncertainty, while the set for the variance 
is deemed as the wind power variability. This approach differs 
from the uncertainty sets on a random variable in 10 by intro¬ 
ducing uncertainty sets on distribution parameters describing 
the random variable. We derive these sets with a data-driven 
approach such that the resulting RCC-OPF remains tractable 
B, B, yielding a practical approach for modeling WPG. 
Although the RCC-OPF model is nominally more conservative 
than the CC-OPF model, its solution incurs a lower operating 
cost when tested against actual realizations of WPG H . 


2 ) Wind Speed Variability: For any subinterval r of hour- 
long interval t wind speed Wr can be written as Wr = 
na, where pt and e^- are the hourly-average wind speed and 
intra-hour, zero-mean variability around p^t Q, respectively. 
To parametrize Wr, we calculate the average wind speed (pt) 
for each hourly interval within the studied period using the 5- 
minute resolution Goodnoe data 0. Next, we bin the hourly 
intervals based on their hourly-averaged wind speed pf This 
binning is in the range from 0 m/s to 25 m/s (the typical 
operating range of wind turbines) and has the width of 1 m/s. 
For each bin we obtain a histogram that represents an empirical 
intra-hour wind speed distribution f^{wr\pt) conditioned by 
the hourly-average wind speed (pt). Each histogram is then 
fit to the Gaussian distribution, yielding estimated mean (p*) 
and standard deviation values (a*) for every bin. Both the 
histogram and its Gaussian distribution are shown in Fig. [T^) 
for several values of pt. Figure [T})) shows that a* scales 
linearly with //* as cr*(/i*) = 0.231+0.197-/i* for p* G [0, 25] 
m/s, which is consistent with velocity distributions for high 
Reynolds number atmospheric fiows 0. If there is no error 
in a given hourly-averaged wind speed forecast p{, the wind 
speed variability is estimated using the linear mapping cr'^ip'^) 
shown in Fig. [Th). The resulting distribution is then given 
as N[wr;p{, (cr*(/i{))^], and can be used in 0 and 10 for 
uncertainty sets accounting for wind speed variability. 


3) Wind Speed Uncertainty and Variability: Wind forecast 
errors cause pt ^ p{. The wind forecast error ct is calculated 
from the NOAA 0 data as et{AT) = p^ — p( (AT), where 
p( (AT) is the forecast for hour t made AT hours in advance. 
The empirical distribution f^(epAT) for AT = 1 hour is 
shown in Fig. [2^). We propose to represent f^{') with a 
generalized Gaussian distribution, /^(•)- 


f^{euAT,p-,p+) = ^ 


dpN[et;p,(T^{p)] 


( 1 ) 


where cr*(/i) is the fit from Fig. [Tp) and the dependence of 
/i+ and p~ on AT is suppressed. The advantage of f^{') 
over a single Gaussian distribution is that it improves the 
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II. Methodology 

I) Data: We use historical wind speed measurements at 5- 
minute resolution from the Goodnoe meteorological station in 
the Bonneville Power Authority (BPA) system 0 and one- 
hour resolution wind speed forecasts produced by the NOAA 
Rapid Refresh numerical weather prediction model 0 for 
the same location. The historical measurements and forecasts 
are detrended by using data from the same calendar season 
(December-February) and range of hours (00:00-04:00 AM). 



Fig. 1: a) Empirical distribution /^(ict|/i*) (histograms) and their 
Gaussian best fits for a few p*. b) Empirical relationship between 
CT* and p* and its linear fit cr*(/i*). 
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a) Wind speed forecast error distribution (e^) 
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Fig. 2: a) Empirical distribution f^{et; AT) (histogram), the best fit 
single Gaussian distribution (black) and the best fit distribution f^{-) 
(red) from Eq. ([T]) for AT= 1 h. b) it{~ and versus AT for 
/t( = 10 m/s. c) cr*(/i{“) and cr*(/i{^) for the data in b). 



multiple turbines by using an aggregated wind power curve 
ma. The conversion of the range [/i{ , is given by: 

^ (3) 

Figure [SJ d) shows and corresponding to p{~ 

and respectively, from Fig. [Jb). In Fig. [2j)), the growth 
of p{p{~^) at larger AT is eventually clipped by p{w) as 
enters Region III of the turbine curve. If the entire range 
is in Region III, then the range [p{p{~),p{p{~^)] 
collapses to zero width around the maximum output. 

The conversion of the range [cr*(/i{~), cr*(/i{^)] is shaped 
by the RCC-OPF formulation in (H. Note that ap can 
be obtained yp{/J^t) ^ using the relation 

cr*(/i) from Fig. [TJ)) and the slope s of the turbine curve 
as cFp = s{iit) ' cr*(/it). However, using the wind turbine 
power curve from Fig. [3^) results in difficult nonconvexity in 
distributionally robust formulations (e.g., Eq. (27)-(28) in BTl). 
To avoid this nonconvexity, we assume that the range on the 
mean value in (Eq. [3]) and range on the standard deviation are 
independent. Therefore, the range on the standard deviations 
is given by: 

(j- = min s(M)cr*(M) (4) 

(jj" = max s(p)(7*(/i) (5) 

Eigure [3p) shows the conversion to the robust interval on the 
standard deviation of WPG for the data in Pig. ^). 



Prediction time (AT), h Prediction time (AT), h 


Fig. 3: a) Typical wind power turbine curve p{/t) with four nominal 
operating regions lO. b) Ranges for the on the mean WPG (i.e. 
forecast error) computed using Eq. for the data in Eig. |2t). 
c) Ranges on WPG standard deviation (variability) computed using 
Eqs. (ID) and ^ for the data in b). 

goodness-of-fit to empirical data lH by encapsulating a linear 


III. Conclusion 

We have presented a data-driven method to develop robust 
intervals for distribution parameters, which preserves the phys¬ 
ical relationship between instantaneous and hour-average wind 
speed and power. This method is suitable for uncertainty sets in 
the RUC © and RCC-OPF O, g|, which can be interpreted 
in terms variability and uncertainty of WPG. The case study 
in i) shows that the RCC-OPF model with these uncertainty 
sets outperforms several benchmarks in terms of several cost 
and reliability metrics. 
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by solving the optimization problem: 

, + \ . a^m2 


arg mm 

fi- 


IAT, 11-,^+) - AT)\^det, (2) 


which minimizes the mean square difference between /^(•) 
and f^{') and ensures a better fit to the historical data than a 
single Gaussian distribution, as shown in Fig. [2^). We interpret 
: range where — 

4 + /i+, .... —....... 

wind speed. Fig. displays 

for different AT. Using cF''{iJi') from Fig. [it), we compute 
the bounds on a* as cr*(/ij~) and cr*(/i{^), which are shown 
in Fig.[T]:). The ranges and [cr*(/i{~), cr*(/i{^)], if 

converted to wind generation as explained below, can be used 
in the RCC-OPF formulation from O, 01. 


the range where and = 

p{ + /i+, as the bounds of the uncertainty set for the mean 

and p{~^ for p{ = 10 m/s 


4) Conversion to Wind Power: We illustrate the conversion 
to wind power using the single wind turbine power curve p{p) 
shown in Fig. [3^), lO. This procedure can be generalized to 
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